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Abstract 

The data of direct single-photon measurements of the WA98 cohaboration in 
the reaction Pb(158 AGeV) + Pb are analyzed within a thermal model with 
a minimum number of parameters adjusted to the dilepton data obtained by 
the CERES and NA50 collaborations in the reactions Pb(158 AGeV) + Au, 
Pb. The agreement of our model with the WA98 data points to a unique large 
thermal source emitting electromagnetic radiation observable in both the real 
and virtual photon channels. 
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I. INTRODUCTION 



Direct real and virtual photons are now widely considered as very important penetrating 
probes of the early, hot and dense stage of relativistic heavy-ion collisions (cf. [0] for a recent 
survey). Efforts have been devoted recently to understand the dielectron spectra 

measured by the CERES collaboration in the reaction Pb(158 AGeV) + Au in the low 
invariant mass region and the dimuon spectra measured by the NA50 collaboration in the 
reaction Pb(158 AGeV) + Pb in the intermediate mass region [Q. Both experiments point 
to a significant excess of dileptons above the conventional sources such as electromagnetic 
decays of light hadrons (resulting in the hadronic cocktail), and correlated semileptonic 
decays of open charm mesons, and the Drell-Yan process. It is commonly believed that, 
in particular, the CERES data, obtained with lead and sulfur projectiles, could indicate 
a drastic change of the properties of hadrons, mainly the p meson, in a hot and dense 
surrounding hadronic medium For instance, transport simulations including in- 
medium modifications of the light vector meson masses and widths are found to account for 
the excess of dielectrons. At the same time "traditional" hydrodynamical calculations 0, 
even when supplemented by the in-medium effects, give a dilepton yield which appears to be 
below the experimental data in the most important invariant mass region, 0.4 GeV < M < 
0.6 GeV. Only if one employs chemical off-equilibrium effects, e.g., parameterized by a finite 
pion chemical potential [|l^, the agreement of the hydrodynamical model calculations with 
data can be improved f^. 

In a recent analysis we proposed a model with a minimum set of parameters, which 
describes at the same time the CERES and NA50 dilepton data by a thermal source. In 
spite of its schematic character, the model is based on the nontrivial observation, made 
recently as a result of rather involved evaluations of the properties of hadrons near the 
chiral symmetry restoration, that the dilepton emission rate from a hadron gas at given 
temperature is fairly well described by the quark-anti-quark annihilation rate at the same 
temperature in a wide range of invariant masses extending to small values M . This is called 
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a " duality" of the hadronic and quark-gluon degrees of freedom in describing the emissivity 
of strongly interacting matter. Such a simplification makes it possible to parameterize the 
thermal source by a Boltzmann like exponential function with effective temperature T^g and 
a normalization factor A^^efr which reflects the space-time volume occupied by the thermal 
source [Q. The most remarkable result of such a parameterization is that the dilepton 
spectra in both the low-mass region, as measured by the CERES collaboration, and the 
intermediate-mass region, as accessible to the NA50 collaboration, can be well described by 
the two unique parameters Tes and Nes- Also the transverse momentum spectra are well 
described . This clearly indicates a common thermal source seen in different phase space 
regions and different dilepton channels. 

In the present paper we employ the above mentioned model for the thermal source to 
analyze the direct real photons in the reaction Pb(158 AGeV) -|- Pb. In line with the hadron- 
quark "duality" we assume that deconfined and hadron matter shine equally bright in a 
temperature region centered around the critical temperature of chiral symmetry restoration. 

An important information, which can be extracted from the photon spectra, is related 
to the transverse collective flow of matter. Unfortunately, the corresponding transverse 
momentum {Q±) spectra of dileptons do not give so far such an opportunity, since in the 
large Q± region, where the flow effect is strongest, the total yield is dominated by the 
background contributions. 

The primary aim of the present note is to show that a thermal source, with parameters 
adjusted to the CERES and NA50 dilepton data, results in a photon spectrum which, when 



including the background, is in agreement with the data of the WA98 collaboration |[TT|] for 
the reaction Pb(158 AGeV) + Pb. 

II. MODEL 



For the emission of photons from a thermal source we employ the rate obtained in |T2 
for the free quark-gluon plasma at temperature T 
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dN 5 2 f P'^ 



1 + ; 



(1) 



asT 

where p'^ = {E,p) = {p±cosh.y,p±sinh.y,p±) denotes the four-momentum of the photon 
with transverse momentum p± and rapidity y; u'^ is the four-velocity of the medium; the 
constant k is according to k, = 2.912/(47r). We use units with h = c = 1. The rate 
(|I|) takes into account the lowest order QCD processes, i.e. the qq annihilation and the 
Compton like process. For the strong coupling strength we take as = 0.3. As shown in 
111, at T = 200 MeV this rate is very close to that of hadron matter when including a 
quite complete list of reactions of light mesons and the decays of vector mesons. This gives 
an additional hint for extending the mentioned hadron-quark "duality" to the real photon 
sector, at least for temperatures in the region of the chiral symmetry restoration. (At much 
higher temperature, the rate in [|13| is probably more appropriate; cf. also [p!4| .) 

Since the thermal photon spectra are thought to be sensitive to the medium's flow we 
consider here the simplest case of a spherically symmetric expansion with velocity profile 
v{r) = VQr/R{t), where r stands for the radial coordinate and R{t) is the radius of the source, 
so that for a constant velocity parameter vo the size R{t) increases linearly with time t. In 
line with our model of the thermal dilepton source we replace the temperature by an 
effective average temperature Tes, being constant. Within such an approximation one can 
factorize the space-time volume V4 and the radiation emissivity. Performing the space-time 
integration in eq. (|l|) one gets the photon spectrum as 

dN 

E— = V,F{E,%s,vo), (2) 

3 /• 3^ dN 
4^J d^x{d^p/E)' 

where s = r/R(t) and V4 = ^ / dtR(tY. Below we use the uniquely fixed parameters 
Teff = 170 MeV and V4 = N^s = 3.3 x 10*^ fm'^, which have been previously adjusted to 
the CERES and NA50 dilepton data. 
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III. ANALYSIS OF PHOTON SPECTRA 



We focus on the direct photon production assuming that the background, related to 



decays of secondary hadrons, is removed from the corresponding data To describe 

the spectra one also needs the contribution from hard initial processes, like the Drell-Yan 
process in dilepton production, which is expected to dominate in the high-p^ region. The 
hard photon yield is generated in our study by the event generator PYTHIA ||T^ with 
structure functions MRS D-', default cut-off parameter p±min = 1 GeV for photons, and 



intrinsic transverse parton momentum spread y {k\) = 0.8 GeV, adjusted to the transverse 
momentum spectra of dileptons in the Drell-Yan region in otherwise the default switches 



are used. Higher order calculations are presented in ||16| , but we try to absorb them into an 
appropriate K factor. 

To check the reliability of estimates of the hard component, we compare in fig. 1 the 
PYTHIA results with the photon data obtained in pp collisions by the E704 collaboration 
at ^ = 19.4 GeV (left part) and by the UA6 collaboration at y/s = 24.3 GeV p| (right 
part). To reproduce the data one has to include the K factors of 3.2 (1.8) for ^/s = 19.4 GeV 
( y/s = 24.3 GeV). As seen in fig. 1, PYTHIA describes then well the data, i.e. the slopes 
of the spectra, which are in the high-p^ region, p_i_ > 2.5 GeV. This conclusion is important 
for the subsequent analysis of the photon data in heavy-ion collisions, where the absolute 
normalization of the hard production (i.e., the number of nucleon-nucleon collisions) is a 
priori unknown. 

The comparison of our thermal model yield plus the hard yield from PYTHIA with the 



WA98 data |Tl[ is displayed in fig. 2 for the case of neglecting the flow, i.e., vq = 0, and in 
fig. 3 for the flow parameter vq = 0.3. (Similar to the PYTHIA calculations we simulate the 
thermal yield eqs. (||, ^ by a Monte Carlo procedure and apply the detector acceptance.) 
As seen in fig. 2, the hard photon production process dominates in the high-p^ region, where 
we have adjusted the PYTHIA calculation to data. At p± < 2.5 GeV one can observe the 
increase of the experimental yield above the hard yield. In this region secondary processes 
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become important. One has to stress, however, that in the region of not too large values 
of p± any hard photon production calculation based on perturbative QCD is not longer 
reliable. In particular, in the region of smaller p± the PYTHIA results become sensitive to 
the cut-off parameter p±min- 

As can be expected from eq. (0), the contribution from the thermal source strongly 
depends on the transverse flow. The moderate value of vq = 0.3, used in fig. 3, is consistent 
with the analysis of transverse hadron flow at kinetic decoupling ||19|. For such a flow 
parameter one finds good agreement with the data fill] when adding the thermal yield and 
the hard contribution (see fig. 3). 

To make more firm conclusions on the role of thermal photons in Pb + Pb collisions 
at CERN-SPS energies, one needs a more reliable procedure to fix the hard background. 
Similarly to dileptons, an accurate adjustment of the hard rate at the high-p^ fail requires 
an improvement of the data statistics. Notice that our present up-scale factor from PYTHIA 
simulations of pp collisions to heavy-ion data can be considered as an upper bound for the 
hard radiation. Nevertheless, a remarkable space is left for the secondary radiation. 

An additional insight can be gained by an analysis of non-central collisions which should 
allow a link to pA collisions, where also data are at hand PDI. This will be subject of a 
separate future study. 



IV. SUMMARY AND DISCUSSION 

In summary we analyze the direct photon production in the reaction Pb(158 AGeV) 
-|- Pb by using a model with a minimum parameter set, i.e., the effective temperature 
and the space-time volume of the thermal source, both ones adjusted to dilepton data in 
similar central reactions. The model employs the hadron-quark "duality" for the rate of 
electromagnetic radiation off matter. We have shown that our model, supplemented by the 
hard photon yield as described by PYTHIA, is in good agreement with the WA98 data. This 
result supports the assumption of an extended and long-living source of the electromagnetic 
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radiation, which can be seen in both the real and virtual photon channels. The comparatively 
large value of the space-time factor V4 in our model can be related to previous expectations 
on a long-living fireball or a phase space overpopulation of hadrons P,p!0[|. 

It should be emphasized that our effective temperature parameter Tes = 170 MeV is 
in perfect agreement with the temperature parameter needed to describe hadron species 
ratios |^. Since Tes is to be considered as average of the temperature, one concludes that 
the electromagnetic probes indeed point to temperatures above the expected deconfinement 
temperature. This corroborates the expectation of exciting deconfined matter in central 
heavy-ion collisions already at CERN-SPS energies. 

The photon spectra are shown to be useful in extracting information on transverse flow. 
For firm conclusions, however, the hard photon production processes must be reliably acces- 
sible. Our approach can be contrasted with attempts to interpret the data either without 
any hard contribution or by the hard yield alone by tuning parameters. Transport 

approaches should smoothly interpolate between these extreme cases. 

We expect that the future analysis of the starting experiments at the relativistic heavy- 
ion collider RHIC at Brookhaven National Laboratory should deliver a higher value of T^s 
since the estimated maximum temperatures will be significantly larger. 
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FIGURES 




r "I 

2 2.5 3 3.5 4 4.5 5 5.5 6 
Pt [GeV] 

FIG. 1. Comparison of PYTHIA results with direct photon data in pp collisions. Left part: 
^/s = 19.4 GeV and E704 data (full dots) 0' "ght part: y/s = 24.3 GeV and UA6 data (full 
squares) [ll8| . 
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FIG. 2. Comparison of the thermal photon yield (thin curve, no transverse expansion) and the 



hard background (dashed curve) with the direct photon data |11| for the reaction Pb(158 AGeV) 
+ Pb. 
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FIG. 3. As in fig. 2, however with transverse flow {vq = 0.3). The solid line depicts the sum 
of the thermal and hard yields. 
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